Signatures of Planck-scale interactions in the cosmic microwave background? 
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Based on a rather general low-energy effective action (interacting quantum fields in classical curved 
space-times), we calculate potential signatures of new physics (such as quantum gravity) at ultra- 
high energies (presumably the Planck scale) in the anisotropics of the cosmic microwave background. 
These Planck-scale interactions create non-Gaussian contributions, where special emphasis is laid 
on the three-point function as the most promising observable, which also allows the discrimination 
between models violating and those obeying Lorentz invariance. 
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Introduction In our present standard model of cos- 
mology (see, e.g., 0), the anisotropies observed in 
the cosmic microwave background (e.sr by the recent 
WMAP and the future PLANCK ^ mission) orig- 
inate from quantum fluctuations of the inflaton field, 
which were amplified and stretched by the cosmic ex- 
pansion during inflation (a very early epoch in the evo- 
lution of our universe). The expansion of our universe 
acts as a cosmic microscope, i.e., these observed modes 
can be traced back to extremely short length scales 
and thus ultra-high energies (such as the Planck energy 
Mpi « 10 19 GeV) which are completely inaccessible with 
present or future terrestrial experiments (e.g., particle 
accelerators). This observation entices the question of 
whether one might observe signatures of new interactions 
(such as quantum gravity) beyond the standard model of 
particle physics in the anisotropies. 

Several different approaches for such expected devia- 
tions from our known laws of physics at highest ener- 
gies have been investigated already. One possibility to 
change the physics at small distances is a modification of 
the dispersion relation (of the inflaton field) implying a 
breakdown of the Lorentz invariance (which would then 
be just a low-energy symmetry), see, e.g., It was 
also shown that an explicit short-distance cutoff has a 
potentially measurable effect on the spectrum Possi- 
ble back-reaction effects resulting from UV-modifications 
have been considered in and a possible change of the 
propagator for the inflaton while preserving the local 
Lorentz invariance has been studied in 0. Several au- 
thors have investigated different effects that result from 
a modification of the commutation relations at high ener- 
gies such as = ihSij(l+rjp 2 ), see, e.g., Finally, 
the impact of different choices of the initial vacuum state 
have been considered in [ij], for example. 

These references already illustrate the importance 
of considering the correct initial state. (One could obtain 
basically arbitrary results by changing the initial state 
correspondingly.) Furthermore, the aforementioned in- 
vestigations were mainly devoted to linear fields, whereas 
one would expect new physics such as quantum gravity 
to impose nonlinear corrections, i.e., interactions. There- 



fore, we shall consider a rather general nonlinear low- 
energy effective action for the inflaton field (interacting 
quantum fields in classical curved space-times) and cal- 
culate the corrections due to new physics at ultra-high 
energies in the anisotropies of the cosmic microwave back- 
ground. For simplicity, we assume that the Planck scale 
is the characteristic scale for new physics from now on 
(but our results can easily be generalized). 

The Model During slow-roll inflation, the space-time 
can be approximated by the de Sitter metric (h = c = 1) 
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with the proper time r and the conformal time t, respec- 
tively, and the Hubble constant f) governing the dynam- 
ics of the scale parameter a(t). Assuming locality and co- 
variance (non-covariant terms will be discussed later) , the 
most general low-energy effective action for a scalar field 
can be represented in a gradient expansion and reads 

C = A(<p) + (dip) 2 B(ip) + (dip) 4 C(ip) + 0[{dpf] , (2) 

with (dip) 2 = (d fl ip)(d fl p) and arbitrary coefficient 
functions A(ip), B{ip), and C(<p), Assuming stability 
B(ip) > and adopting the usual field re-definition 
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The first two terms correspond to the usual scalar field 
theory whereas the third contribution of this low-energy 
effective action is a remnant of new physics (e.g., the in- 
teraction with massive gravitons) at ultra- high energies. 
Identifying the scalar field <f> with the inflaton field, we 
assume the potential for chaotic inflation V(4>) — m 2 (j> 2 , 
but our main results can easily be generalized to other po- 
tentials (see discussion below) . Now we split the inflaton 
field <f> into a classical background field 4>q plus its quan- 
tum fluctuations <j>(t,r) = 4>o(t) + Q(t,r). In the usual 
slow- roll approximation m <C -f), the slowly decaying clas- 
sical background solution reads 4>o{t) oc exp(— m 2- r/3.f)), 
and the Lagrangian for the fluctuations is given by 
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with g = 2W{(j) )dcf> /dT » -20 O W (</> )m 2 / '(3£) denot- 
ing a roughly constant effective coupling strength. Note 
that the interaction term (<9</>) 4 in the Lagrangian © also 
implies small corrections such as <E> 2 and thus a small re- 
definition of the propagation speed c e g, the mass m, and 
the amplitude of the inflaton field <E> etc. As usual during 
inflation, the mass term to 2 $ 2 can be neglected for the 
fluctuations <!> due to m <c Sj (terms like A<E> 3 will be 
discussed below). Since we are mainly interested in the 
three-point function, we omitted the higher-order terms 
C[$ 4 ], but the associated four-point contribution can be 
calculated via the same formalism (see discussion below) . 
For the sake of brevity, we omit stating the higher-order 
corrections 0[(dcj)) 6 } + 0[$ 4 ] + 0[g 2 ] from now on. 

Hamiltonian After a Legendre transformation, we ob- 
tain the low-energy effective Hamiltonian density 



H 



n 2 



a 2 (V$) 2 



9 



n 3 



n 



(v$) ; 



(5) 



The cxplicitc time-dependence of this Hamiltonian can 
partly be removed via the (time-dependent) canonical 
transformation n — > an and <&—><!>/ a generated by the 
squeezing operator S, giving rise to an additional term 
to the Hamiltonian: In the squeezed representation with 
the Schrodinger equation id(S \^))/dt = H S (S |*)), the 
Hamilton density reads 

n s = \ [n 2 + (v$) 2 ] - ^ [*n + n$] 

+4 [n 3 - n(v$) 2 ] = Ho + n n + n g . (6) 

This form allows for an intuitive physical interpretation: 
Due to the cosmic expansion a — e^ T , the modes are 
being continuously stretched and hence the effects of the 
interactions H g oc 1/a 2 decrease with time whereas the 
influence of the expansion Hp, oc a increases. 

Initial state Since the interactions are strong at early 
times, the derivation of the correct initial state requires 
the treatment of strongly interacting quantum fields in 
curved space-times, which is a nontrivial task, cf. [llj| . 
Let us consider the evolution of the modes in more detail: 
Due to spatial homogeneity, the co-moving wavenumber 
k corresponding to a spatial dependence of e lk v remains 
constant - but the physical wavenumber k/a decreases 
(stretching). Hence modes with different k leave the 
Planck scale at different times, i.e., when k/a = 0[MpJ. 
Assuming that the rate of the cosmic expansion given 
by the Hubble parameter .fj is much slower than the 
internal evolution rate of these Planckian modes (scale 
separation fj <C Mp\), the most natural initial state is 
the adiabatic vacuum state, cf. During the further 
expansion, the wavelength increases until k/a <C Mp\ 
holds and the low-energy effective action becomes 
valid 111 view of the adiabatic theorem, the modes 
evolve adiabatically and hence stay in the adiabatic vac- 
uum as long as the external time-dependence due to the 



cosmic expansion is slow compared to the internal dy- 
namics. If we now choose an intermediate time to such 
that $j <C k/a{to) <C Mpi, the effects of both, the cosmic 
expansion Hfy and the interaction H g are small and can 
be treated as perturbations 01 As a result, the quan- 
tum state (adiabatic vacuum) of these modes at that time 
to can be derived via the adiabatic expansion [liij 
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where \^ n ) are the instantaneous eigenvectors of the 
Hamiltonian © with appropriate phases and No — 
1 + 0[,f> 2 ] is a normalization. For g = 0, the Hamiltonian 
Ho + Hf, can be diagonalized exactly (i.e., to all orders in 
Sj) via a Bogoliubov transformation Cfc = a^bf, + /JfcoLfc- 
Up to second order in Sj (the accuracy we are inter- 
ested in), the Bogoliubov coefficients otk and ftk read 
a k = 1 + iaS)/2k and /3 k = iaf)/2k - a 2 ^ 2 /Ak 2 . The 
remaining corrections can be calculated with stationary 
perturbation theory 
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where (H + H%) = E° |*°) arc the unperturbed 
eigenstates obtained via the aforementioned Bogoliubov 
transformation. 

Final state The impact of the cosmic expansion H^ 
on the modes increases with time and their evolution fi- 
nally becomes non-adiabatic when k/a — 0[Sj], i.e., when 
they cross the horizon and freeze. Therefore, in order to 
calculate the final state, we switch from the Schrodinger 
representation (thereby undoing the canonical transfor- 
mation S etc.) to the interaction picture: The Heisen- 
berg operators carry full dynamics for the unperturbed 
problem Hq + H^ to all orders in 9) (cf. 0) 
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Note that a k annihilates the adiabatic vacuum for 
vanishing g, which can be used as a consistency check. 
The remaining small interaction Hamiltonian H g acts on 
the quantum state, which enables us to calculate the final 
state via time-dependent perturbation theory 



|out(t)> = | X—* y dt'H g (t) ) |in(io))+C% 2 ]. (10) 

to 

As another test of the consistency of the used formalism, 
one may check that the resulting final state |out(t)) is 
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indeed independent of to as it should be. Since we are 
mainly interested on modes which freeze well before the 
end of inflation, we may approximately extend the inte- 
gration to infinite proper time r = oo which corresponds 
to vanishing conformal time t = 0. 



Three-point function With the final state obtained 
after performing the integration, one can calculate the 
three-point function of the frozen fluctuations and ob- 
tains (after some algebra) the following spectrum 
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with the abbreviations k + = k\ + k 2 + k 3 and ej = kj/kj. 
The remaining permutations (denoted by "perm.", i.e., 
fei — » fe, — > ^3, and fc 3 — > fci etc.) have to be included 
as well. The spectrum is scale invariant as one would 
expect since the sequence: leaving of Planck scale and 
stretching — > horizon crossing and freezing occurs on all 
scales (just at different times). The first contribution 
3/fci. on the second line results from the n 3 -term whereas 
the remaining ones appear because of the n(V$) 2 -terni 
in the interaction Hamiltonian H g . 

An estimate of the relative size of this non-Gaussian 
contribution can be obtained by comparison to the 
(Gaussian) two-point function 

^H^= = 0[gSj 2 } = 0{foW(Mm 2 f)} « 1 • (12) 

Note that m can be eliminated with the classical equa- 
tions of motion and the Friedman equation, i.e., gfi 2 = 
0{(j> W((f>o)m 2 f)} = 0[W(<p )^ 3 My<j)o]. Usually the 
typical scale of new physics such as quantum gravity is 
expected to be the Planck scale. Therefore, one would 
naturally assume that the coupling behaves according to 
W ~ 1/Mpj if W did not depend on </> . In this case, 
the relative size of the effect would be suppressed with 
i^ 3 /M|j and thus probably hard to measure. 

However, one should also bear in mind that the usual 
picture of inflation involves a large value of the inflaton 
field 4>o S> Mpi and hence W((j)o) could be much bigger 
than 1/Mpj. Similarly, the scale for new physics might 
be below the Planck scale or the inflaton potential V(<j>) 
could differ from m 2 <j> 2 . In all these cases, the effect 
could be far stronger than fj 3 /Mpj. On the other hand, 
the employed formalism only works if gf) 2 <C 1 holds, i.e., 
if the relative size of the effect (|12|l is small compared to 
one (which is consistent with the observations |2|). 

Generalizations Let us discuss the impact of 
momentum-independent interaction terms such as A$ 3 
[in contrast to contributions like (9<i>) 4 ] which occur if 
the potential V is not just given by the mass term 
V(4>) 7^ m 2 (j> 2 . In contrast to H g , these momentum- 



independent interaction terms become important at late 
times. Hence the three-point spectrum generated by the 
A<f> 3 -term can be derived from time-dependent perturba- 
tion theory and behaves as fc+ / (fcifc 2 A: 3 ) 3 . This spectrum 
is dominated by large length scales and thus not scale 
invariant - which should make it easy to distinguish be- 
tween the effects of the momentum-independent interac- 
tion A<& 3 and those generated by high-energy (Planckian) 
interactions such as the spectrum (|1 1|) . 

The derivation of the four-point function can be ac- 
complished by means of the same formalism - but ne- 
cessitates higher orders in the adiabatic expansion etc. 
The relative size of the non-Gaussian contribution (com- 
pared to the Gaussian part of the four-point function, 
for example) can be estimated to 0[W(0o).f> 4 ]- Hence 
the comparison of the relative sizes of the non-Gaussian 
contributions of the three-point and the four-point func- 
tions, respectively, yields the unknown ratio Sj(f)o/Mpy 

Lorentz invariance So far, we assumed that the 
Planckian interactions respect the local Lorentz invari- 
ance and hence started from a covariant action J2J . How- 
ever, it may well be that the usual local Lorentz invari- 
ance is just a low-energy symmetry and broken at the 
Planck scale. In this case, many more terms can occur 
in the low-energy effective action - even if we still de- 
mand spatial homogeneity and isotropy according to the 
cosmological principle: First of all, bilinear terms such 
as $ 2 and (V 2 <I>) 2 exactly correspond to a change in the 
dispersion relation mentioned in the introduction. As- 
suming adiabaticity (which may not by valid for all dis- 
persion relations), we may treat those corrections with 
the presented formalism and it turns out that they ba- 
sically generate a global shift of the two-point spectrum 
to lowest order, which is probably hard to measure. 

The next terms occurring in an expansion into powers 
of $ and the numbers of derivatives are exactly the terms 
n 3 and n(V<i>) 2 discussed before jj^. However, without 
Lorentz invariance, the pre-factors of these two terms are 
no longer necessarily equal. Therefore, measuring and 
comparing the relative strength of the two contributions 
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in the above three-point spectrum (fTTfl provides a pos- 
sibility to detect violations of the Lorentz invariance at 
the Planck scale. 

Summary Assuming locality and covariance, we de- 
rived the most general low-energy effective action of the 
scalar inflaton field (including possible remnants from 
Planckian interactions) via a gradient expansion. Based 
on this ansatz, a formalism for the treatment of this 
(strongly) interacting quantum field theory in a curved 
space-time and the derivation of the quantum state via 
the adiabaticity assumption was presented. This enabled 
us to predict the induced spectrum of the three-point 
function - which is quite robust against uncertainties in 
cosmic evolutions and can clearly be distinguished from 
other effects such as a A$ 3 -term - making it a promis- 
ing observable. This concrete prediction of the explicitc 
spectrum facilitates a refined search in the present/future 
WMAP and PLANCK data with good statistics and 
hence high sensitivity. Abandoning covariance, it might 
even be possible to detect violations of the Lorentz in- 
variance at the Planck scale. 
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